Nonsymmetrized Correlations in Quantum Noninvasive Measurements 
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A long-standing problem in quantum mesoscopic physics is how to order operators in quantum 
noise expressions like (J(— where is the measured current at frequency uj. Symmetrized 
ordering of operators describes a classical measurement while nonsymmetrized ordering corresponds 
to a quantum detector, sensitive to either emission or absorption of photons. We show that both 
ordering schemes can be embedded in quantum weak-measurement theory taking into account a 
finite interaction time between the system and the detector. Importantly, the nonsymmetrized 
ordering reveals its nonclassical nature already in second-order correlations, contrary to symmetrized 
ordering. This feature can be related to the squeezing of the many-body state of the transported 
electrons in a ac-driven tunnel junction. 



Quantum measurement theory has been based on the 
projection postulate [T], although nowadays it includes 
generalized schemes based on auxiliary detectors, de- 
scribed mathematically by positive operator-valued mea- 
sures (POVM) [2 . To specify a POVM requires argu- 
ments based on physical considerations such as detector 
efficiency, or the assumption of thermal equilibrium. A 
real physical interaction generally leads to backaction on 
the system to be measured, which makes the interpre- 
tation of measurements difficult. Hence, all detection 
schemes are in general invasive as the measured system 
is perturbed. The disturbance is strongest for projective 
measurements, as the information in the measurement 
basis is completely erased. In contrast, other POVM 
schemes can be much less disturbing, as is often the case 
in experiments [3H5]. 

To avoid invasiveness, Aharonov, Albert and Vaidman 
[Hj studied the limit of weak measurement, in which the 
system is coupled so weakly to the detector that it re- 
mains almost untouched. The price to pay is a large de- 
tection noise, which is however completely independent 
of the system. The gain is that other measurements on 
a non-compatible observable can be performed. After 
the subtraction of the detector noise, the statistics of the 
measurements has a well-defined limit for vanishing cou- 
pling, which for incompatible observables turns out to be 
described by a quasiprobability and not a real probability 
distribution [TUB]- 

The most common weak-measurement theories assume 
that the system-detector interaction is instantaneous 
[SHU]- Such a Markovian measurement scheme justi- 
fies the symmetrized ordering of operators, relevant for 
many experiments [3]: ~> v(~ w)J(o;) + 

I(u))I(—U)))/2 (denoting quantum expectation values 
(X) = TiXp for an initial state p). Here, I(w) — 
J dt I(t)e lujt is the Fourier transform of the time- 
dependent current I(t) in the Heisenberg picture. How- 
ever, certain experiments, e. g. involving absorptive de- 



tectors, are well described by nonsymmetrized correla- 
tors like (I(-uj)I(uj)), corresponding to emission noise 
for uj > [U 15 -20 . These experiments clearly lie be- 
yond the scope of Markovian weak measurement theory. 

In this Letter, we formulate a general theory of weak 
detection which allows for the description of nonsym- 
metrized correlators. Whereas no Markovian scheme of 
weak measurement predicts emission noise, we show that 
the latter appears naturally if one allows for a measure- 
ment with memory. In fact, non-Mar kovian weak mea- 
surements follow just from a few natural assumptions 
imposed on the POVM in the limit of weak coupling. 
By further requiring that no information transfer oc- 
curs in thermal equilibrium we fix uniquely the scheme, 
which contains the detector temperature as only param- 
eter. Our scheme interpolates between emission and ab- 
sorption measurements, when varying the detector tem- 
perature. As it is independent of other properties of 
the detector, our scheme applies to circuit QED, meso- 
scopic current measurements and quantum optical sys- 
tems equally well. Interestingly, applied to a simple har- 
monic oscillator, the scheme reveals classical equations of 
motion but with the initial distribution corresponding to 
the Glauber-Sudarshan P-function [3T]. Contrary to the 
instantaneous measurement, the non-Markovian scheme 
can violate weak positivity |19j . To test it, we propose 
a measurement of photon-assisted current-fluctuations, 
which are shown to violate a Cauchy-Schwarz type in- 
equality. Identifying the finite-frequency current oper- 
ators with quadratures in analogy to quantum optics, 
we show that the thus created non-equilibrium state of 
the current is squeezed and therefore has essentially non- 
classical correlations. 

We start by developing a the general framework of 
weak quantum measurement based upon the POVM for- 
malism including non-Markovian features. We consider 
a set of n independent detectors continuously recording 
n time-dependent signals dj{t) for j — l,...,n. Each 
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detector is related to an observable Aj. For example, 
n ammeters are inserted in a complex circuit: cij(t) is 
the recorded current in the branch j and Ij(t) the cur- 
rent operator in that branch. Note that in general Aj (t) 
and Ak{t') do not commute even if j ^ k since Aj and 
Ak may not commute with the hamiltonian. We want 
to relate classical correlators of measured quantities like 
(ai(t\) ■ ■ ■ a n (t n )) to their equivalent for weak quantum 
measurements (■ ■ - ) w - These should involve linear corre- 
lators of the Aj , which can be taken at different times to 
allow for memory effects of the detectors, while preserv- 
ing causality. The requirements of linearity and causal- 
ity are fulfilled by replacing aj(t) in the correlator by a 
superoperator J dt'A* - * (t') and perform time ordering, 
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Here T denotes time ordering with respect to the argu- 
ments in brackets, p is the density matrix, and A are 
superoperators defined as: 
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t—t' 



(f ) = 9j (t t')A%t') + fj(t - t')A%t')/2 . (2) 



The superoperators A 



c/q 



act on any operator X like 



an anticommutator/commutator: A c jX — {Aj,X}/2 and 

AjX — i\Aj,X\. In the above expressions we supposed 
for simplicity that the detectors are in a stationary state 
so that only time differences tj — t'j matter. It is easy 
to show that the correlator ([!]) is invariant if we change 
the order of the aj , and that all combinations of Aj that 
preserve causality can be reproduced. 

We will also assume that the average of single mea- 
surements coincides with the usual average for projec- 
tive measurements, i.e. (aj(t)) w — (Aj(t)). This implies 
cjj(t — t r ) = S(t — t'). Other choices of g simply mimic the 
effect of classical frequency filters. Thus the only freedom 
left is the choice of the real function fj that multiplies 
Aj. Note that fj{t) can be non-zero for t > without vi- 
olating causality, since it is accompanied by Aj and only 
future measurements are affected. For the last measure- 
ment, future effects disappear because the leftmost A q 
vanishes under the trace in Eq. ([I]). For simplicity, we 
will assume a single f = fj, independent of j. The limit 
/ = corresponds to the Markovian case. 

Now we want to show that correlations obeying these 
requirements can be obtained from the general quantum 
measurement formalism. Based on Kraus operators K 
[23] , the probability distribution of the measurement re- 
sults is p = (K) for KX = KXK^ , where the only condi- 
tion on K is that the outcome probability is normalized 
regardless of the input state p. Here we need K to be 
time-dependent. In general, we assume that .?C[j4,a] is a 



functional of the whole time history of observables A{t) 
and outcomes a(t). We shall assume that the functional 
K is stationary so it depends only on relative time argu- 
ments. 

The essential step to satisfy Eq. ([I]) is to take the limit 
K r*> 1 which corresponds to a noninvasive measurement. 
This can be obtained from an arbitrary initial POVM by 
rescaling K[A, a] — > K v = C(ri)K [r]A, rja] with rj — » 0, 
which defines p n = (K v ). Here C(rj) is a normalization 
factor. 

The desired correlation function ([T]) can be derived 
by the special limiting procedure for an almost general 
POVM, namely 



(oi(ti) • • • a n (t n ))t 



lim(ai(ii) ■ • • a n (t n )) v (3) 

Tf— >-0 



where the average on the right-hand side is with re- 
spect to p v . We assume the absence of internal cor- 
relations between different detectors, namely K[A, a] = 
l~Y[j K[Aj,aj], where T applies to the time arguments 

of A. 

Expanding 



K[A,a]/k[a] = 1 + / At'F[a, t']A{t') + <D(A 2 



(4) 



wc find 



K/\k[af 



At' (2ReF[a,t']A c (t') 



-ImFla^'jA^t 1 )) +0(A 2 ) 



(5) 



Here, k[a] is a functional amplitude independent of the 
properties of the system which represents the detec- 
tion noise. As we want the measurement to be non- 
invasive to lowest order, we impose the condition that 
/ F[a, t'] \k[a] | 2 Da vanishes; Da is the functional measure. 
Our conditions imply that J 2a(t)ReF[a,t']\k[a]\ 2 ~Da = 
S(t-f), and we get f(t-f) = J 2a(t)luiF[a 1 t'}\k{a\\ 2 T>a. 
Thus, the most general weak Kraus operator takes the 
form given in Eq. ([5]), which is our main result. A partic- 
ular Gaussian example of a POVM realizing this scheme 
is presented in the Supplementary Material A. 

To discuss the consequences of different forms of /, we 
now calculate the noise spectral density, 



S ab {w) = / dt e iut (a{t)b(0)) t 



(6) 



An important special case is a system in a thermal equi- 
librium state, p ~ exp(—H/kBT). We further assume 
that the averages of A and B vanish. If the detec- 
tor temperature Td is equal to T and in the absence 
of other nonequilibrium effects (like a bias voltage, or 
special initial conditions), we expect that no informa- 
tion transfer from the system to the detector occurs, i.e., 
that S a b(Lu) = 0. This requirement leads to a necessary 
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condition on the form of / (see Supplementary Mate- 
rial B): f(w) = -i{2n B {u) - 1) = -% coth(huj/2k B T d ), 
where n B (uj) is the Bose distribution at temperature T d . 
Equivalcntly, f(t) = — k B T& coth(irtk B Td/h)/h (at zero 
temperature f(t) = —l/irt). We use the name equilib- 
rium ordering for this special choice of /. In a recent 
work, the zero temperature case has been called time- 
normal (24] . 

The necessary form of / is also sufficient. Indeed, 
the property S a b(u) = follows from the fluctuation- 
dissipation theorem [25] 

df e iut (A(t)B(0)) T = J dt e lult+huj ' kBT (B(Q)A(t)) T , 

(7) 

because for an arbitrary stationary state we get 

S ab (u) = J dt e lut {e hui / 2kBTd B{Q)A{t) 
-e~ hul/2kBTd A(t)B{0))/smh(hu}/2k B T d ) , (8) 
which reduces to the emission noise [18 



'dt (0(-u)A(t)B(O) + 0(u)B(p)A(t)) (9) 



for zero detector temperature Td — 0. Thus, for Td T, 
S a b(ijj) is in general not equal to zero and contains infor- 
mation about the system. Reversing the sign of / will 
reverse U) in ([9|, transforming it into absorption noise. 
Equilibrium ordering is relevant for experimental situa- 
tions like in [3] and consistent with the quantum tape 
[T§] or photodetection model [2U] if the temperature of 
the tape (or the photons) is T d . 

It is interesting to note that for this special choice 
of / the higher-order fluctuations also vanish if p oc 
exp(— H/k B T) and T = T d . We can write the Fourier 
transform of ([TJ as 



Wfe] 



4_„±hu k /2k B T A±(j. \ 



fc ± 



2sinh(^w fe /2A: B T) 



with A+X = AX and A~X = JtA. Now, we can split 
p = p 1 / 2 /) 1 / 2 , expand the above expression as a sum of 
operator products and move one factor p 1 / 2 leftwards and 
the other rightwards so that they meet again at the trace 
sign, which gives (JlJ in the form 



J d n t e^^^TrpTjJ 



(11) 



E 



±e ±hu */ 2kBT A%(t k T ih/2k B T) 
2smh{hu k /2k B T) 



Shifting t -> t ± ih/2k B T and using TiA q ... = leads 
to 

J d n t e^*' lUktk TrpT]jA q k (t k )/2ismh(hu; k /2k B T) = 0. 

k 

(12) 

The proof is hence completed for all nonzero frequencies. 
In every physical situation it extends also to zero fre- 
quency because the cumulants are continuous functions 
of frequency. 

Surprisingly, the equilibrium ordering differs qual- 
itatively from the symmetrized one, when one con- 
siders weak positivity, i.e., the classical nature of 
second-order correlations [SI [26]. The symmetrized 
correlation matrix C a b = (ab) = (AB + BA) /2 
is positive definite so the Gaussian probability dis- 
tribution oc exp(— ^ af) C~ f) 1 a6/2) reproduces all first 
and second order symmetrized quantum correlations. 
This is not the case in equilibrium ordering. If 
T d > 0, then already for a system in equilibrium 
at T = and u > we get (a(-uj)a(u)) = 
~(A(u)A(-u))e- hu / 2kBTd /2smh(huj/2k B T), which is 
negative. For Td = and stationary states again the 
weak positivity holds, because the correlation matrix ^ 
is positive definite. However, it is violated for nonstation- 
ary states. This can be demonstrated using a two level 
system with the Hamiltonian H = hfla z /2, with observ- 
able A = a x + a z and the initial state p(0) = (1 + & y )/2. 
By direct calculation (Supplementary Material C) we find 
(a 2 (0)) lu = — (2/tt) lnfiioo, where t^ is a cutoff set by in- 
trinsic decoherence or detector backaction. Thus, weak 
positivity is obviously violated. In reality (a 2 (0)) w corre- 
sponds to (a(O)fo(O)) where the detectors a and b monitor 
the same observable A = B in ([lj. 

There is an interesting connection between equilibrium 
ordering and the Glauber-Sudarshan P function [21] . Let 
us take the harmonic oscillator H = hVl (p 2 + x 2 ) /2, 
with [x,p] — i, and consider correlations with respect 
to the quasiprobability ([TJ) with antisymmetric f(t). In 
this case the time ordering is irrelevant as shown in Sup- 
plementary Material D. Let us define A(t) = A c (t) + 
/ dt' f(t - t')A q (t')/2. The evolution of x(t) and p(t) is 
just classical, 



x(t) = xcos(Ot) + psu\(VLt), 
p(t) = pcos(Q,t) — xsm(Vlt) 



(13) 



while for t = (since f(t) is antisymmetric and real) we 
have 



x = x c + ip q f(n)/2, 

p = p c -ixif(n)/2. 



(14) 



We now define the ladder operators through a) — (x 
ip)/y/2 with [a, of] = 1. This leads to a = a c - f(n)a q /2 
and a) 



/(f2)a' 9 /2. In the zero-temperature case, 



and defining a 



ip)/y/2 we get the 
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quasiprobabilistic average (a n a* k ) = Tra n pa) k . On the 
other hand, this is a property of the Glauber-Sudarshan 
function P(a), defined by p = J d 2 a P(a)\a) (a\ for nor- 
malized coherent states a\a) — a\a), (a\a) = 1 [2"T] . 
Since (a n a* k ) P = J d 2 a a n a* k P(a) = Tra n ptf k , we find 
that the initial quasiprobability for a zero-temperature 
detector is identical to P(a,t — 0). It is interesting to 
note that reversing the sign of / leads to the Husimi- 
Kano Q function instead of P [2], while / = gives the 
Wigner function [71127]. 

The fact that we obtain the P-function shows the deep 
connection between the non-Markovian weak measure- 
ment formalism and the quantum-optical detector the- 
ory. One of the interesting consequences is that zero- 
temperature equilibrium ordering is consistent with pho- 
toabsorptive detection schemes, in which the P-function 
appears naturally [2J. It is also interesting to draw a 
link between the violation of weak positivity in equilib- 
rium ordering and the properties of squeezed states. The 
ground state of a harmonic oscillator fulfills (x 2 ) = 1/2, 
which corresponds to (a; 2 ) p = 0. A squeezed state can be 
such that (i 2 ) < 1/2, still minimizing the Heisenberg un- 
certainty principle. This translates into a negative vari- 
ance of the position described by the (quasiprobability) 
P-function, i.e. (x 2 )p < [28] and is therefore equivalent 
to a violation of weak positivity. 

Let us now consider how our results apply to the case of 
current fluctuations in mesoscopic conductors The quan- 
tum description of the noise in the junction, Si(lo) = 
J dt e iut (SI(t)SI(0)), where 61 (t) = I{t) - (/(<)), will 
depend on the choice of / in 0. For / = 0, we get 
symmetrized noise Sf = Gh^2 ± w(lu ± eV/h)/2, where 
G is the conductance, V is the constant bias voltage and 
w(a) — acothfia/kpT) [TJ]]. For a choice of / with an 
arbitrary T<j, we obtain Si — Sf — Ghwd{u)- Hence, the 
detection schemes differ by a term that is independent of 
the voltage and the temperature of the system, making 
it impossible to detect non-classicality in this scheme. 

An experimentally feasible test of squeezing and vi- 
olation of weak positivity is possible using a coherent 
conductor (e.g. a tunnel junction for the sake of sim- 



plicity) subject to an AC voltage bias V(t) 
Consider the classical inequality 



Vac cos Vlt. 



\8I{u)-5I{-uj)\ 2 > (\SI(lu)\ 2 ) > Re(5I 2 (uj)) . (15) 
For symmetrized ordering one gets [29] 

{{SI(u), SI(uj')}) /2 = 2irhGj2 s (" + w ' - 2mn ) 

m 

^ J n (eV ac /KSi)J n -2 m (eV ac /ft£L)w{u - nil) , (16) 

n 

where J n are the Bessel functions. In the case of equi- 
librium ordering at = one only has to subtract 
2tt1iG\ui\6(ui + ui') from the above result. As shown in 
Fig. [l] the classical inequality is violated for uj = VL in 



Re(/ 2 (fl)) 

</(-n)/(n)> 
<{/(-n),/(n)}>/2 




eV a 



3.o m 



FIG. 1: (color online) The plot shows the quantum correlation 
functions in tunnel junctions normalized by 2nGhQ,to at T = 
0. The emission nois e, at Td = 0, (red line) violates the 
classical inequality ( 15 1 for a certain range of eKc/ffl (shaded 
region) . The violation is equivalent to the squeezing condition 
(JTrj) for the symmetrized noise (blue line). 



a certain range of eV ac /Ml, but only in equilibrium or- 
dering. This can be reinterpreted in terms of the exis- 
tence of squeezing in the quantum shot noise : consider 
the two quadratures associated with the finite-frequency 
current operator: A = i\8I[uS) — 5I(—ui)]/2 and B = 
[5I{uj)+5I{~ui)]/2. Using {[I{u)),I{-u))}) = 2t Ghuj 1 we 
find [3D] ([A, B]) = itoGhu (with the total detection time 
to). Thus the squeezing condition [2] 



(A 2 )<\([A,B})\/2 



(17) 



is related to the violation of weak positivity, (A 2 ) w < 
in equilibrium ordering ([I]) with Td = and allows to 
violate Eq. (15). Hence, according to Fig. [I] quantum 
shot noise with ac-driving creates current states, which 
resemble squeezed light for a certain range of the ac- 
voltage. 

In conclusion, we have presented a theory of a generic 
weak-measurement scheme that includes emission noise. 
It requires a non-Markovian POVM with a specially cho- 
sen memory function /, which has no analog in the 
Markovian picture. The scheme is consistent with the ex- 
pectation about the absence of information flow between 
system and detector in equilibrium at a given tempera- 
ture. Hence any detection requires a nonequilibrium sit- 
uation. Another direct consequence is that even the sim- 
ple Markovian detection process must involve a nonequi- 
librium detector state. Finally, nonsymmctrized order- 
ing leads to a violation of weak positivity, which can be 
tested experimentally by violation of suitable inequali- 
ties, equivalent to the squeezing condition in some cases. 
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Supplementary Material 
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A. GAUSSIAN NON-MARKOVIAN POVM 



An example of a POVM leading to a non-Markovian weak measurement is the Gaussian Kraus functional 



K[A, a] oc Texp / dt 



-(A(t) - a{t)) 2 + / dt' if(t - 0(2a(t) - A(jt)0(t - t'))A(t') 



(A.l) 

Here, the first term in the exponent is the Markovian part, while the second term describes the non-Markovian mea- 
surement process including a fixed but arbitrary real function fit), characterizing the memory effect. The Heaviside 
function 6 is necessary to ensure the normalization of the Kraus operator. Following the standard procedure we find 
the Kraus superoperator in the form 



K[A, a] oc 
Texp 



(A.2) 



2(A c (t) - a(t) f + (i 9 (<)) 2 /2 + / dt' 2f(t - t')(a(t)A q {t') - 6{t - t')(A c (t)A q {t') + A q {t)A c (t')) 



dt 

To prove the normalization, jDa(K) = 1, we perform the Gaussian integral over a (time ordering is no problem if 
kept up throughout the calculation) and get 



dt(A q {t)) 2 /2 + / dt'6(t' - t)2f(t - t')A q {t')A c {t)dt' 



Da K = T exp 



-0(t - t')2f{t - t')A q {t)A c (t')dt' + f{t - t')f(t - t")A q (t')A q (t")dt'dt"/2] , 



(A.3) 



where we have ordered properly A q it') and A c it). In the power expansion, omitting the identity term, the leftmost 
superoperator is always A q . Since TiA q ■ ■ ■ = we obtain J Da (K[A, a]) = 1 or J Da K^K = 1. In general, we define 



K[A, a] for n measurements as if [A, a] = Tn,'-S"[^ 



Note that putting A — gives Gaussian white noise p oc 



e 2a ~ , which leads to the large detection noise in the weak limit, p v oc e 2ri a , that has to be subtractcd/deconvoluted 
from the experimental data. 



B. FIXING THE MEMORY FUNCTION / 



Since the detector function /(w) should be system-independent in thermal equilibrium, any system can be used to 
determine it. We therefore consider a 2-level system with A = B = a x and H = Mla z /2. The requirement S(uj) = 
is equivalent to 



Re 



e«*dt ((1 - if(uj))a x (t)a x (0) + (1 + ifiuj))a x (0)a x it)) = . 



The equilibrium state reads p = (1 — <j z t&n\i(hil/2kBT))/2 and 

(v x (0)vx{t)) = (v x i-t)a x iO)) = cos(Qt) +it&nh(hn/2k B T)sm(nt) , 



and (JB . 1 [) leads to the requirement that 
Re ' 



I 



e + iu + ifl e + iuj — iVt 



1 



1 



e + iuj + ?0 



- m 



tanh 



2k B T 



(B.l) 



(B.2) 



(B.3) 



vanishes for e — > + . Since = + nS(x), we can ignore the delta function for u ^ ±f2, and the vanishing of 
(B.3) reduces to Re/(cj) = . As / must be independent of the system, we are free to choose f2, Re/(w) = must 
hold for all ui, including u> — ±57. So / is purely imaginary, and (B.3) reads 



S(uj + 6(uj - ft) - (5(u + ft) - 5(lu - ft)) tanh 
Therefore, 1 ± Im/(w) tanh(ftft/2yfc B T) at u = ±ft, and f(u) 



( hfl 



\2k B T 
icothifuv/2k B T) 



W(w) 



(B.4) 
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C. VIOLATION OF WEAK POSITIVITY 

From wc find for f(t) = -l/nt 

(a(t)b(s)) w = ({A(t),B(s)}/2 + f idt'[A(if),B(s)]/2n(t-tf)+ f i&s'[B{s'),A(t))/2itts-s>)) . (C.l) 

\ J — oo J — oo / 

For H = Mla z /2 and A = B = a x + a z we find A(t) — a x cosQ,t — a y sin fit + <t z and i[^4, A(i)]/2 = <7.j.smf2i — 
a z sin fit + cfy(cosQt — 1). Therefore, for p(0) = (1 + <J y )/2 we get 

(a 2 ) =2+-/ dt . (C.2) 



For small t the integral is convergent but for large t only cos flt/t converges. The remaining integral J dt/t diverges 
logarithmically and one should put a cutoff at t^. 

Certainly no experiment will record infinite correlations. The cutoff is in practice bounded by the decoherence 
time of the system and the measurement noise (which also diverges). The infinity would occur only in the limit of 
zero measurement strength and a perfect two-level system, which is impossible. 

D. HARMONIC OSCILLATOR 

For the harmonic oscillator with H = p 2 /2m + moj 2 x 2 and [x,p] = ih, we have [x(t),p(t')} = ihlcos(il(t — £')), 
[x(t),x(t')] = — ihl sin(f2(i — t'))/mfl, [p(t),p(t')} = —ihls'm(il(t — t'))mf2, so the commutator depends only on the 
difference t — i', which applies also to superoperators. To see that the time ordering is irrelevant, let us take linear 
functions A,B of x and p and calculate 



A c (t) + J dt' f(t-t')A q (t')/2j [B c {s)+ J ds' f(s- s')B q (s')/2 

T (A c (t) + f dt' f(t - 0^(0/2 ] (b c (s) + f ds' f(s - s')Bi(s') ) (D.l) 



= / dt' f{t-t')6{ s -t')[B c {s) 1 A q {t')]/2 + J ds' f(s - s')6(t - s')[B q (t),A c (s')}/2 
= ( duO{u){f{u + t-s)[B c {s),A q {s-u)]+ f{s-t-u)[B q (t + u) 1 A c {t)])/2. 



The last expression vanishes because of the antisymmetry of / and the fact that the commutators depend only on the 
difference in time arguments. The proof generalize to multiple products because [B c ^ q (s), A c / q {u)] is proportional to 
identity superoperator for A,B — x,p. 



